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GLOBAL SOLUTIONS FOR THE ZERO-ENERGY 
NOVIKOV-VESELOV EQUATION BY INVERSE SCATTERING 


MICHAEL MUSIC AND PETER A. PERRY 


Abstract. Using the inverse scattering method, we construct global solutions 
to the Novikov-Veselov equation for real-valued decaying initial data qo with 
the property that the associated Schrodinger operator —d x d x + qo is nonneg¬ 
ative. Such initial data are either critical (an arbitrarily small perturbation of 
the potential makes the operator nonpositive) or subcritical (sufficiently small 
perturbations of the potential preserve non-negativity of the operator). Previ¬ 
ously, Lassas, Mueller, Siltanen and Stahel proved global existence for critical 
potentials, also called potentials of conductivity type. We extend their results 
to include the much larger class of subcritical potentials. We show that the 
subcritical potentials form an open set and that the critical potentials form 
the nowhere dense boundary of this open set. Our analysis draws on previous 
work of the first author and on ideas of P. G. Grinevich and S. V. Manakov. 


Contents 


1. Introduction 

2. Preliminaries 

2.1. The 9-Problem 

2.2. Properties of the scattering transform 

3. Topology of the set of subcritical potentials 

4. Continuity of Maps 

4.1. Continuous dependence of t on q 

4.2. Continuous dependence of reconstructed q on t 

5. Symmetries of the scattering transform 

6. Solving the Novikov-Veselov equation 

7. Conclusion 
References 


1 

6 

6 

7 

Jjj 

1C 

1C 

12 

10 

19 

21 

21 


1. Introduction 

The Novikov-Veselov equation at zero energy is a (2 +1 )-dimensional completely 
integrable system given by 

d T q = 3 x q + d x q - 3 d x (uq) - 3 d x (uq) 

( 1 . 1 ) 

d x u = d x q, 
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where 

i/dd\ _1 (d . d\ 
x 2 V&zq +l dx 2 ) ’ x 2 VS X1 l dx 2 )' 

This equation is part of a family of equations, parameterized by the energy E, 
studied by A. P. Veselov and S. P. Novikov !28) because of its connection with 
isospectral flows for the Schrodinger operator at fixed energy E (in our equation, 
E = 0). In their original work, Veselov and Novikov considered periodic solutions; 
P. G. Grinevich, S. V. Manakov, R. G. Novikov, and S. P. Novikov developed the 
inverse scattering method for nonzero energy E and initial data vanishing at infinity 
[sans mini. We refer the reader to Kazeykina [14] for a review of results and 
results on long-time asymptotics at nonzero energy, and to [6] for a comprehensive 
review of the literature and for a self-contained introduction to the inverse scattering 
method for the NV equation at zero energy. 

To study two-dimensional inverse scattering problems, Faddeev [8] introduced 
the complex geometric optics (CGO) solutions. In our case, these are a family of 
solutions to the zero-energy Schrodinger equation with spatial asymptotic behavior 
parameterized by a complex parameter k. Throughout this paper we will treat 
x = (xi,x 2 ) £ K 2 as a complex number (aq + ix 2 ) when multiplying so that 
kx = (k\ + ik 2 ){x\ + ix 2 ). Assuming q G L P (R 2 ) for a fixed p G (1, 2), we say that 
ip(x, k) is a CGO solution for k G C \ {0} if 

- d x d x ij) + q(x)ip = 0 

e~ ikx ip(x,k) - 1 G W l ’ p { R 2 ) 
for some p G (1,2), where p is the Sobolev conjugate: 

1 _ 1 1 

P P 2' 

It is more convenient to work with the normalized solutions 

p(x, k) = exp(—ikx)ip(x, k ) 

which satisfy 

(1.2) d x (d x + ik)p(x , k) = q(x)p(x, k ) 
where 

n{ - ,k) — 1 G W /1, ^(R 2 ). 

If, for some fc, equation d admits a solution h(x) with h(x) G IV 1 ’P(IR 2 ), unique¬ 
ness fails and k is called an exceptional point. Below we will describe a spectral 
condition on q(x) which rules out exceptional points. 

The associated scattering data t for the potential q is 

(1.3) t(fc) = / ek{x)q(x)p(x, k) dm{x) 

J R 2 

where dm( ■) denotes Lebesgue measure on R 2 and 

ek(x) = exp i (kx + kx). 

The map T : g —> t defined by m and m is the direct scattering map. 
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The normalized CGO solution p.(x, k) can be recovered through the scattering 
data, and the potential can in turn be recovered from the reconstructed solutions 
/r(x, fc). Indeed, if 


s( k) = 

7 rk 

the function /i(x, k) satisfies the 9-equation 


(1.4) 


(dfc/z) (x,k) = e- x (k)s(k) n(x, k) 

Kx, -)-ieL f (€) 


where r £ (jp, 2) and r is the Sobolev conjugate of r. The potential q(x) is recovered 
in turn using the reconstruction formula 

(1.5) q(x) = — d x ( / e- x (k)s(k) fj,(x,k)dm(k) 

K \J C 

The map Q : t —>• q defined by m and (USD is the inverse scattering map. 

Formal computations [3] (see also 10 §3.6]) show that, if q{x,r) solves NV, then 
the function t (k, r) = T(q (-, r))(fc) should evolve by the linear evolution 

(1.6) t(fc, r) = exp | 'ir (j; 3 + ^ t(fc, 0). 

Putting these ingredients together-the direct scattering transform T, the lin¬ 
earization m, and the inverse transform Q-we get the solution formula 


(1.7) 


q(x, t) = Q 


ex P(*((' ) 3 + (-) 3 t) (Tq 0 )(-) 


(x)- 


From dEl and it follows that, to solve the NV equation, we need to solve 

the 9-problem with parameters 


( 1 . 8 ) 

and recover 


(d k p) (x, k, t) = e- x (k)s(k, r) /x(x, k, r) 
p(x, •, r) — 1 e L r (C) 


(1.9) q(x, t) = % -d x yj e- x (k)s(k 1 r)/x(x, k, r) dm(k)j 

where we define 

n ^ t ( fc > r ) 
irk 

As we will show, global existence for the formal solution ( 11 . 71 ) within a class of 
real-valued, regular, and decaying initial data is guaranteed by a simple spectral 
condition on the initial data. For a real-valued function q € Lf oc (R 2 ) for somep > 1, 
we say that the operator L = —d x d x + q is nonnegative if, for all <p £ C° 


( 1 . 10 ) 


\(d x ip)(x)\ + q{x) |<p(x)| dm(x) > 0. 


Murata [T8J defined a trichotomy of potentials for the zero-energy Schrodinger 
operators (see Simon [25] for an earlier and closely related trichotomy in the study 
of Schrodinger semigroups). This trichotomy plays a crucial role in our analysis 
and may distinguish distinct dynamical regimes for the NV equation. 


Definition 1.1. Fix p > 1. A potential q £ Lf oc (R 2 ) is 
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(i) subcritical if there is a positive Green’s function for the operator L, 

(ii) critical if there is no positive Green’s function but L > 0, and 

(iii) supercritical if —d x d x + q is not nonnegative. 

Murata (see m Definition 2.1] and accompanying discussion) shows these are 
the only three possibilities, so the condition ( 11 . 101 ) implies that q is either critical or 
subcritical. Theorem 2.12 of [P7] implies that if q is uniformly locally in L P (R 2 ) for 
p > 1, then — A + q > 0 if and only if there exists a positive distributional solution 
to the Schrodinger equation (we use here the fact that, for p > 1, the uniformly 
locally L p potentials are contained in Kato class J \ 2 considered in [7j; see [TJ §1.2, 
Remark 1(b)] and see [7J, §2.5] for further discussion and references to the original 
work of Allegretto, Moss, and Piepenbrink). Additionally, if q(x) has sufficient 
decay, the equation —d x d x ip + qip = 0 admits a positive, bounded, distributional 
solution %p in the critical case, and a positive solution with logarithmic growth in the 
subcritical case m Theorem 5.6]. In Section [3] we show that subcritical potentials 
are an open set in the topology L p ( R 2 ) with p £ (1,2) and p > 1 + 2/p and the 
critical potentials are the boundary of this set. 

The asymptotic behavior as |x| —> oo of the positive solution is closely related, 
via a Fourier-like duality, to the asymptotic behavior as |fc| 4- 0 of the scattering 
transform t(fc). Indeed, Nachman [21] Theorem 3.3] showed that, for critical poten¬ 
tials, |t(fc)| < C\k\ e for some e > 0 and small k. This means that s(k) = O (l^l 6 ” 1 ) 
is in Lf oc (C) for any p £ (l,oo) and the integral operator (12.81) for the problem 
ed has range in L°° functions using (12.31) of Lemma l2Jl By contrast, Music [19] 
showed that, for subcritical potentials, 

s(fc) ~ —=---fl + O ((— log(|fc| 2 ) 1 ') N ) 

V |k|->o k log(|/c| 2 ) V l V \ J )) 

and the same key integral operator does not have range in bounded functions. This 
singularity is the cause of some technical pain in our computations but, thanks 
to a density argument (see Lemma 931 and Corollary 15.31) . does not pose essential 
problems. In Proposition 15.41 we prove that any scattering transform that has this 
singularity and also has certain decay properties (see Definition 12.71) must come 
from a critical or subcritial potential. 

To state the theorem, denote by L P (R 2 ) the space 

{/ £ L p (R 2 ) : (1 + | • | 2 ) P/2 /(■) € L p (R 2 )} 

and by IF™’ P (R 2 ) the space of measurable functions with n weak derivatives in 
L p (R 2 ). 

Theorem 1.2. Let qo £ WA p (K 2 ) forp £ (1, 2) and p > 1 be a critical or subcritical 
potential. Then q(x, r) given by (11.71) is a global-in-time, classical solution of NV 
(ED with lim T _»o t) = qa(x) pointwise. 

The regularity assumption guarantees that the inverse scattering method will 
yield a classical solution of NV. By contrast, Angelopoulos [1] recently proved local 
well-posedness for NV in the space LP(R 2 ) for any s > 5 . The inverse scattering 
method captures global properties of the solution at the expense of more stringent 
regularity and decay assumptions on the initial data. 
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The condition that the potential is subcritical or critical is necessary in the above 
theorem using our methods. This is because in |20j . Music, Perry, and Siltanen con¬ 
struct radial, compactly supported, supercritical potentials qo where the scattering 
transform t(fc) has a circle of singularities and the current formulation of the in¬ 
verse scattering method fails. P.G. Grinevich and R.G. Novikov m give explicit 
examples of point potentials with similar contour-type singularities. Taimanov and 
Tsarev (see |27j and references therein) construct examples of supercritical poten¬ 
tials that give rise to solutions of the zero-energy NV equation that blow up in 
finite time. 

The technical core of our result is a careful analysis of the 9-problem (11.81) with 
parameters x, t. The first technical challenge is that the scattering data t (k) for 
our class of potentials may be mildly singular as |fc| l 0 (see Lemma l2i>l) . We 
deal with the singularity by proving continuity properties of the direct and inverse 
scattering maps and approximating t(k), the scattering transform of a potential q 
in our class by a smooth function with compact support away from k = 0. 

The second technical challenge is that, in order to prove that the reconstructed 
potential CD2) solves the NV equation, we must first prove that it is a real-valued 
function. We will show that this is the case by adapting ideas of Grinevich [9] and 
Grinevich-Novikov m- 

Lassas, Mueller, Siltanen, and Stahel m study the same problem for initial 
data of conductivity type which have very well-behaved scattering transforms as 
shown by Nachman in his 1996 paper [21: on Calderon’s problem. These authors 
obtain a number smoothness and decay results for this class, but they cannot prove 
that the inverse scattering method yields solutions to the Novikov-Veselov equa¬ 
tion. Combining our result with the earlier one from Lassas, Mueller, and Siltanen 
m gives global solutions to Novikov-Veselov which satisfy the decay condition 
\q(x,r) | < ( x)~ 2 if we start with a critical potential q(x, 0) £ C^R 2 ). Although 
this special case has already been proved by Miura map methods by Perry [23| , the 
proof give here is arguably more direct and simpler. 

In Section O we state some theorems that will be used throughout the paper. 
In the first half, Section [2711 we recall facts about the d problem, and in Section 
12.21 we recall results of Music [19] on the properties of the scattering transform for 
critical and subcritical potentials [T9]. Music’s characterization of the range of T 
will allow us to analyze the inverse map Q rigorously. In Section [3j we prove that 
the set of subcritical potentials is open and that the set of critical potentials is its 
nowhere dense boundary. In Section (T[ we establish continuity in q for the direct 
scattering map T, and continuity in s(fc) of the inverse scattering map Q. In Section 
El we exploit ideas of P.G. Grinevich and S.V. Manakov [10] to prove necessary and 
sufficient conditions for the scattering data t (k) to be the scattering data of a real 
potential q(x) repairing a gap [Ml in the proof of [TB] Theorem 4.1], We apply 
the ideas of Grinevich-Manakov m directly to the zero-energy case instead of 
treating it as a limit of negative-energy inverse scattering (see also Grinevich 0)- 
With this result, we prove that the reconstructed potential q(x, r) stays real for all 
time. Crucially, the reality of the inverse map is needed to show that the evolved 
/j,(x,k,T) continue to solve the problem (11.21) . In Sectional we follow the method 
outlined in the review article of Croke, Mueller, Music, Perry, Siltanen and Stahel 
[6] to show that the potential solves the Novikov-Veselov equation. We discuss some 
open problems in section [3 
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2. Preliminaries 

2.1. The 9-Problem. Here we review some important facts about the d operator 
and its inverse that will be used in what follows. Further details and references 
can be found, for example, in [22j section 2] or |2| chapter 4]. We denote by P the 
integral operator 

i p f] ( z ) = - [ —- 7/(0 <MC) 

71 Jc z ~ C 

where, in the integrand, z = z± + iz? and C = Ci + *( 2 - We recall without proof 
the following basic estimates. In what follows, C p (resp. C Pj q) denotes a numerical 
constant depending only on p (resp. p, q). 

Lemma 2.1. (i) For any p G (2, 00 ) and f G T 2 p/(p+ 2 )(C), 

(2-D \\Pf\\ p <C p \\f\\ 2p/{p+2 y 

Moreover, VPf G L 2 p/(p+ 2 )(C) and 

(2-2) l|VP/|| 2p/(p+2) <C p ||/|| 2p/(p+2) . 


(ii) For any p , q with l<q<2<p<oo and any f G L p (C) flL 9 (C), the estimate 

(2-3) \\Pf\\oo < C p , q \\f\\ LPnLq 

holds. 

(Hi) If v € L S (C), p G (2, 00 ) and q € (2, 00 ) with q~ x + 1/2 =p _1 + s _1 , then for 
any f G L P (C), 

(2-4) \\P(vf)\\ q <C Ptq \\v\\ s \\f\\ p . 

The operator P is the solution operator for the problem du = f: 

Lemma 2.2. Suppose that p G (2, 00 ), that u G L p ( C), that f G L 2p ^ p+2 \C), and 
that du = f in distribution sense. Then u = Pf. Conversely, if f G L 2p ^ p+2 \ ( C.) 
and u = Pf, then du = f in distribution sense. 


We will need the following generalized Liouville Theorem for quasi-analytic func¬ 
tions, due in this form to Brown and Uhlmann. Music m has an extension of this 
lemma to include certain negatively weighted L p spaces. 

Theorem 2.3 ([5] Corollary 3.11]). Suppose that u G L P (C) PI Lf oc (C) for some 
p G [1, 00 ) and that 

du = au + bu 

for a and b belonging to L 2 ( C). Then u = 0. 


Finally, solutions of du = f for rapidly decaying / have a large-z expansion. 

Lemma 2.4. Suppose that p G (2, 00 ), that u G L p ( C), that f G L 2 n^ p+2 ' > (C), and 
that du = f. Then 




yj +1 


<7/(0 dm{ 0 


j=o 


G L p { C). 


We will also use the following estimates from [4] on the fundamental solution 


_3 

kernel for the linear problem v T = df.v + d x v associated to the NV equation. 
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Lemma 2.5. Let 

( 2 . 5 ) I T {x) = J e ir(k 3 +k 3 )-i(kx+kx) dm ( k y 

Then: 

I t (x)=t- 2 / 3 I 1 (xt- 1 / 3 ) 

and the estimates 

\h(x)\<C(l + \x\)- 1/2 , 

| V x h{x)\ < C 

hold. 

For a proof see [4] Proposition 5.4]. 


2.2. Properties of the scattering transform. The following is a reformulation 
of Theorem 1.4 and Lemma 5.2 from Music m ■ Instead of stating the lemma in 
terms of the space that the potential q lives, we state the assumptions that the 
scattering transform satisfies. This is the purpose of defining the space X^ r in 
Definition 1271 See Remark 12.81 for more details. 

For subcritical potentials, we let ip(x) be the unique positive solution of — d x d x ip+ 
qif = 0 such that 

( 2 . 6 ) ^(x) = alog |x| + 0 ( 1 ). 

We also set 

(2.7) Coo = lim —^ [ (i/>(x) - alog |x|) dm{x). 

r->oo ttK J\ x \ < r 

Note that for critical potentials, we have a = 0 and Cqq = 1, and the following 
results are originally due to Nachman ED Section 3] in this case. 


Lemma 2.6. Let q & L p p ( R 2 ) withpG (1,2), and p > 1, and suppose that q is either 
critical or subcritical. Let if{x) denote a positive function solving (—d x d x +q)ip = 0. 
Then q{x) has no exceptional points, and the scattering transform t(fc) of q{x) has 
the following properties: 

(i) For small k and some e > 0, we have 


t(fc) = 


1 

2 Co o 


7 TCI 

o(log \k\ + 7 ) 


0(|fc| e ), 


where c Q0 , a are defined in equations (12.61) and (12.71) . 

(ii) t (k)/k € L r (\k\ > e) for r £ (p', 00 ) and every e > 0. Moreoever, t (k)/k £ 
C°( C\0). 

(iii) If q £ Wff ,p ( R 2 ) for n> 1 then |fc|"s(fc) £ L r (C). 


In fact, we will use the above properties for scattering transforms coming from 
real potentials to define a class of scattering transforms. 


Definition 2.7. The space Xf r for n > 1, r £ (1,2), and e > 0 is the closure of 
C) functions which satisfy the relation kf(k) = —kf(—k) in the norm 

ll/IU«_ r = ||/||l 2 (C) + IK • )"/( • )llir' + eni r( C ) + ||( • )"/( • )IU r (C)- 
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Remark 2.8. The symmetry requirement guarantees that the reconstructed poten¬ 
tial is real-valued; see Lemmas o and 15.21 of Section [5] The conditions in this 
definition are preserved under the linearized NV flow (tin . By Lemma 12.61 scat¬ 
tering transforms s(fc) coming from potentials q £ Wf’ p (M. 2 ) are in X^ r for all 
r £ ( p ', oo) and e > 0 . 

These properties of s(k) £ X^ r are used in Music T9 in the inverse problem 
to reconstruct /i(x, k) and obtain a large-fc expansion of /x(x,k). For the large-fc 
expansion to exist pointwise, we need the decay k n s(k) £ L 2+e n L 2 ~ e ( C) (see 
Lemmas 5.5 and 5.6 of BSD- We need this decay for all k away from the origin, 
and in this region it suffices to have s £ L 2 (C) instead of L 2+e ( C). This is why we 
only need the homogeneous norm on L r +e (C). With this norm and the embedding 
Xf +1 r C L\, we have / i(x, ■ ) — 1 £ 17(C) and for every x we get the expansion 

(EH). 

In what follows, we define the operator T XtT by 
(2-8) (T x>r f) (A;) = P k [s( •, r)e_ x ( • )J] (, k ) 

where 

[Pkf](k) = - f - -1-/(0dm(C). 

7T J C « - C 

We now recall from m that the solution of (11.81) is differentiable in the param¬ 
eters x and r. 


Lemma 2.9. If s £ X ^ r then the unique solution p,(x, •, r) — 1 £ 17(C) of equation 
CD is a times differentiable in x and m times differentiable in t for n > 3m + |a|. 
Additionally, the derivatives of the map {x,t) /r(x, - ,t) — 1 £ L r (JC) satisfy 

d™Df:iJ,(x, ■ ,t) £ L r (C). The derivatives are given by 

(2.9) d™D«v(x, k, t) = [I — T^-'Pkffx, ■ , T ) 

where 


(2.10) f(x, k, t) = [d™D^,s(k,r)e- X (k)] n(x, k,r), 
and P k [s(k,T)e- x (k)f(x,k,T)\ £ 17(C). 

Finally, it follows from Lemma \2 .41 that /j,(x, k,r) has a large-fc expansion. We 
refer the reader to m for a full proof. 

Lemma 2.10. Ifs£ Xf r and if /i solves the equation (11.81) with ^(x, •) — 1 £ L r (C), 
then /i admits the large-k expansion 

(2.11) v(x,k,T) = l + ^ fflj ^ T ) +o(|fc|-") 

i=i 

for fixed x and r. Moreover, we may take a spatial derivatives and m time deriva¬ 
tives with n > |a| + 3m to get 


d?D°(n(x,k,r)-l) 


n—\ot\—3m 

E 


y 


+ o 



+ |a|+3 



1=1 

We will use the above lemma to take up to four spatial derivatives, or one time 
and one space derivative of ai(ir,T). This requires t(A:) £ Xf r and is why we 
assume that q{x, 0) £ Wp ,p (R 2 ). It is also important to notice that these are not 
weak derivatives as —id 1 q = ai(x,r) and Wp ,p (R 2 ) C C 3 (R 2 ). 
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3. Topology of the set of subcritical potentials 


The set of subcritical potentials is actually bigger than the set of critical po¬ 
tentials. Murata EH Theorem 2.4] shows that the set of critical potentials is not 
stable under compact perturbations. Here we go a step further and show the sub¬ 
critical potentials are an open set in the topology of L P (R 2 ) for certain p and p. 
In the corollary to the theorem we show that subcritical potentials are open in the 
TTp ,p (K 2 ) topology for the same p and p as in Theorem II.21 

Theorem 3.1. The set of subcritical potentials is open in the L P ](R 2 ) topology for 
Pi £ ( 1 , 2 ) and pi > 1 + 2/p and the set of critical potentials is its boundary. 


Proof. We prove that the set of subcritical potentials is open by finding positive 
logarithmically growing solutions, (f v , to (—A + v)<j> v = 0 for all v close enough to 
a given subcritical potential q. 

Let Gq(x) = — 1/(27t) log |x|. Then, by Nachman [2TJ Lemma 3.2] we have for 
/ £ L p j](M 2 ) with P 2 > 1 and p 2 £ ( 1 , 2 ) 


G 0 * f + — log \x\ 


f dm(x) 



and 

l|VG 0 * f\\ L P2 < c||/||lp2. 

Modifying both of these slightly, we let 

Tf = G 0 * / + log(|x| + e) f fdm(x) 


and get 

\\Tf\\ W l,P1Z < c\\f\\ L P2. 

In particular, by Sobolev embedding we have 


\\Tf ||l- < H/ILg. 

Let 4> q be a positive solution to the Schrodinger equation and let c aa be given by 
m - By Nachman [2T] Lemma 3.5] for critical potentials and Music Q15] Theorem 
1.4] for subcritical potentials with Coo 0, we have that [I + Go * (q ■)] is invertible 
for critical/subcritical q on Wjp (M 2 ) for /? > 2/p. If Coo = 0, we may use the same 
argument from Music [T9] Theorem 1.4 and Lemma 4.1] to scale q , and therefore 
(j)q, so that Coo ~f~ 0 and we may reduce to the case. So for some e > 0 the operator 
[I+Go*(v ■)] is invertible for all potentials v satisfying ||g — v\\ l p < e. Thus we have 
distributional solutions to (—A + v)(j> v = 0 by taking cj) v = [I+Gq*(v ■ )] _ 1 Coo. 
Because q is subcritical (j) q = alog|a;| + 0(1) is a positive solution. The solutions 
<p v satisfy the equation 

(fv — ^oo Go * ( V(f)v')' 

Taking the difference between these two solutions we find 

(frq — <t> v = G 0 * (v<f> v - q<f> q ) + log(|x| + e) J(y<j> v - q(f q )dm(x) 

~ ^ lo g(M + e) J[vcj) v - q<j> q )dm{x) 

= T(v<t> v - qcj> q ) - log(|a;| + e) j (' v<j> v - q<j> q )dm(x ) 
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which becomes 

\<t>q ~4>v | < c\\v<f) v - # 9 || L P 2 + log(|a;| + e)||w^ - q</) q \\ L i 

for P 2 £ (1,2) and p 2 > 1- We may choose these constants by taking /3 close enough 
to 2 /p so that p 2 = pi — /3 > 1. Also, we have the embedding ||/||li < c||/|| L P 2 so 

4>v > <t> q - c\\v<t> v - # 9 || l p 2 - ^-log(|x| +e)\\v<t> v -q<j> q \\ L P 2 . 

The right hand side is positive and logarithmically growing for \\v<p v — q(j) q \\ L P 2 
small. The norm is small because \\(f> v — 4> q \\ w i,p < c\\q — v\\ l p by the second 

resolvent formula applied to the operators [I + G 0 * (v ■)]. The function </>„ is 
a positive distributional solution to (—A + v)(f v = 0 so by [7J Theorem 2.12] v is 
critical or subcritical. It follows from Murata [18] Theorem 5.6] that v is subcritical 
because the positive solutions for critical potentials in this weighted space have the 
asymptotics </> = c + o(l) whereas the positive solutions for subcritical potentials 
obey (j) = a log |x| + 0 ( 1 ). 

In Theorem 2.4 of the same paper, Murata showed that nonnegative compact 
perturbations of critical potentials are subcritical and nonpositive perturbations are 
supercritical. Thus critical potentials are a subset of the boundary of subcritical 
potentials. Looking at the form (11.1011 we see that limits of subcritical potentials 
must either be subcritical or critical. Therefore, critical potentials form the entire 
boundary. □ 

Now we apply the above lemma to the set of potentials under consideration in 
Theorem 11.21 

Corollary 3.2. The set of subcritical potentials is open in the W/’ P (M 2 ) topology 
for p £ (1,2) and p > 1 and the set of critical potentials is its boundary. 

Proof. Be Sobolev embedding W/ ,P (M 2 ) C T*(]R 2 ) for all s £ (p,p). In particular 
we may choose s < 2 but close enough to 2 so that p > 1 + 2/s. The result then 
follows from the previous theorem. □ 

4. Continuity of Maps 

In this section, we prove that the scattering data, t, depends continuously on the 
potential q and that the solution of (11.81) is continuous in s (k) in a suitable sense. 
Using the continuity of the inverse scattering map defined by the 0-problem, we 
can show that our reconstructed q{x , t) is differentiable. In fact by approximating 
the scattering transform s (k) by s e (fc) G Cf°( C) in the Xf r norm, we approximate 
the solution q(x,r) by the reconstructed q e (x,r). In this way, the results showing 
that each q e {x,r) solves the Novikov-Veselov equation carry over to the general 
subcritical case so long as the initial data qo has enough derivatives. 

In what follows, B(X) denotes the Banach space of bounded operators from a 
Banach space X to itself with the usual norm. 

4.1. Continuous dependence of t on q. Recall from Nachman EH §i] that if 
p £ (1,2), u £ VU 1 ’P(K 2 ), / £ L P (R 2 ), k £ C \ {0}, and d x (d x + ik)u = / in 
distribution sense, then 

u = g k * f 

for a convolution kernel gk which satisfies the estimate 
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\\9k*f\\ W i,e < co(k,p) H/llp. 

Let us define 

Sk(q)f =gk* (qf )■ 

It follows from the estimates above that for q £ L P (R 2 ), 

(4-1) \\Sk{q)\\ b(w^p) < co (k,p) ||g|| p 

while 5fe(g)l := gk * q satisfies 

(4-2) 11^(9)111^1 ,p < c(k,p) \\q\\ p . 

In (14.11) we used the Sobolev inequality WfW^ < c p H/H^i.p- 
We now consider the continuity of p,(x, k\ q) defined by 

g{x, k;q) -1 = (I - S' fc (g)) _1 5' fc (g)l 

as a function of q. 

Lemma 4.1. Fix p £ (1,2), k £ C\ {0}, and qo £ L P (R 2 ). Suppose that 

R(qo) ■= (I-Skiqo))- 1 

exists as a bounded operator from VL 1,P (K 2 ) to itself. There is a number r > 0 and 
a constant c(p , k, qo) so that for all q £ L P (R 2 ) with \\qo — q\\ p < r, the estimate 

IK'.M) ~ T{ - ,k,q 0 )\\ w x,p < c{p,k,q 0 ) ||g-g 0 || p 

holds. 

Proof. In what follows c denotes a constant depending only on k, p , and qo whose 
value may vary from line to line. We will use the estimates 

(4.3) HS'fc(g) — Sk(qo) llg^up) < c ||q — <?o|| p 

(4-4) 11-5*(9)1 - S'fe(9o)l|| B( ivi.? ) < c ||9 - 9o || p 

which follow from ED , (BUD , and linearity in q. 

By the second resolvent formula, for q sufficiently close to qo, we have 

R(q) - R(q 0 ) = R(q 0 ) ( S(q ) - S(q 0 )) R{q) 

so that 

R(q) = [I- R(q 0 ) ( S(q ) - 5(g 0 ))] _1 R(qo) 

These computations can be justified if ||-R(go) [S{q) — jS(go)] Wb^w 1 ^) < 1/2- From 
estimate ED, we see that it suffices to take 

ll9-9o|| p < (^2c ||R(9o)||g(wi.p)) 

Thus, taking r = (2c ||i?(9o)|| B (v^i,P)) , we have 

(4.5) \\R(q)\\ 

BiW 1 ’*) — 2||-R(9o)|| BiW 1 ’*) ■ 

Using the identity 

p(x, k, q) - p(x, k , q') = R{q)S k {q) 1 - R(qo)S k (q 0 )l, 

the estimates ED , ED , ED, and the second resolvent formula, we recover the 
claimed estimate. □ 


Now we can prove continuity of t as a function of q. 
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Lemma 4.2. Suppose that p £ (1,2), p > 1, that q £ L P (R 2 ), and {q n } is a 
sequence from L P (M. 2 ) withq n —>• q in L p ( R 2 ). Finally, lett n = T(q n ) andt = T{q). 
Then, for all non-exceptional nonzero k, t n (k) —> t(fc) pointwise. 

Proof. Using Lemma |4~T1 we estimate 


|t n (AO - t(fc)| < 


e k {x)(q n (x) - q{x))p n {x) dm(x) 


e k {x)q(x))(p n {x) - p(x)) dm(x) 

< II Qn - q\\Li\\Pn\\L°° + |M|l,i||Mn ~ F\\L°° 

< I \Qn - Q\\l p p (II Fn — lllwLP + 1) + IkllijIlMn “ Mil W 1 -* 
and conclude that t„(fc) —> t(fc). 


□ 


4.2. Continuous dependence of reconstructed q on t. For notational conve¬ 
nience, we combine the time dependence of the scattering data 

s (k, t) = exp(ir(A; 3 + fc 3 ))s(fc, 0) 

with the spatial oscillation e- x (k) into the factor e lrS where 

S(x, k, t) = -h ( fc 3 + fc 3 ). 

T 

We note the following estimates. 

Lemma 4.3. Suppose that s £ L 2 (C) and f > 2. Then ( I — exists as an 

operator in B(L r ) and 


sup 

(x,r)eR 2 xR 


(I-T^)- 1 


B{L f ) 


< 00 . 


Proof. We will prove the assertion in four steps. First, we will show that T X)T is a 
compact operator on L r ( C) with 

( 4 -6) sup \\T x , T \\ B(Ln < C(f) ||s|| L2(c) . 

(i,r)el J xR 


Second, we will prove that ker(/ — T XjT ) is trivial, so that, by the Fredholm 
alternative, the operator (/ — T XjT ) _1 exists for all (x, t) £ R 2 x R. 

Third, we will show that 


(4.7) 


lim sup ||T| t 

T->oo a , eR2j \ t \>T 


B(L f ) 


= 0 


and that for each T > 0, 


(4.8) 


lim sup \\T 2 T \\ =0. 

|t|<T, |x|>R ’ 


It then follows that there are numbers R and T for which || T 2 , 
whenever |x| > R or |r| > T, and that therefore 

, -l 


I B(L f ) 


< 1/2 


(4.9) 




B(L f 


< 2 
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on this set. Using the identity (I — A) 1 = (I — A 2 ) 1 (J + A), the estimates (14. 6|l . 
and (1491) , we recover 


(/ - T x , r ) 


-1 


B(L f ) 


< 2 


(i + c II s IIl 2 (c)) 


for all such (x, r). 

Fourth, we will use a continuity-compactness argument to show that for any 
positive numbers R and T, 


sup 

|x|<R, |t|<T 


(J - r X;T ) _1 


B(L f ) 


is finite. 

(1) The norm estimate is immediate from (12.41) with p = q = r. To show that 
T x>t is compact, it follows from (14.61) and density that we may take s G Cq°(C) 
without loss. It suffices to prove that the Banach space adjoint T' x T = e lrS sP k 
is a compact operator on L r (C). Let f l be the support of s. If / G L r (C) then 
P k f G L 2f ^ f ~ 2 \C) by ([24]) while VP k f G if { C) by (ES). Thus 

\\Pf\\w^'(C) <c(l + |fi| 1/2 ) ||/||-, 

and compactness follows from the Rellich-Kondrakov Theorem. 

(2) To see that ker (/ — T x>t ) is trivial, note that by Lemma [2.21 ip G ker (/ — T XtT ) 
if and only if di/j = e lTS sip. Since s G L 2 (C) and i[) G L r (C), it follows from Theorem 
12.31 that ip = 0. 

(3) We will prove (14.71) and (14.81) by proving the corresponding estimates for 
{T x t ) 2 on L r '(C). First, note that T' x t = e lrS W where W = s P k is a compact 
operator independent of (x,r). We have the estimate 

uniform in ( x , r). For any e > 0 there is a finite-rank operator F on L r (C) so that 
||W — F \< e. The operator F is a finite sum where <pj G L r (C), 

ipj G L r (C), and (•, ■) is the usual dual pairing of L r (C) and L r (C). A density 
argument shows that we may take 'tpj and tpj in 5(C) without loss. Thus, it suffices 
to show that 


(4.10) lirn sup \((p,e iTS i/j)\ = 0 

T ~> 00 i6i 2 . |t|>T 

and that, for each fixed T > 0, 

(4.11) lim sup I (ip, e lTS ip) I = 0 

|t|<T, |x|>R 1 


where ip and tp belong to 5(C). 

Let f(k) = tp(k)ip(k). A short computation shows that, for S = S(x,k,r), 
(p, e lTS ip) = J I T (x - y) (J --1 /) (y) dy 

where I T is given by m and 

(- 77-1 /) (s) = ^ J e- k (x)f(x)dm(k). 

We now appeal to Lemma T2. 5 1 to estimate 

(4.12) | (ip, e irS i/>) | < jT- 2 ^l + \(x-y)/T\)^ 2 \^f)(y)\dm(y). 
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To prove (|4.1Q|) . we estimate the right-hand side of (|4.12|) by T 2 / 3 Wj 7 1 || l i. This 
gives (TO) . To prove (14.111) we make the change of variables £ = {x — y)/ t 1 ^ and 
obtain 


|(^e^}| < /(l + |£ir 1/2 |(J-7) {x — £t 1//3 )| dm(0 
<C(f,N)(h(x,T)+I 2 (x,T)) 


where C(f, N) depends on the Schwartz seminorms of / and 

h(x,r)= f (l + ^|)- 1 / 2 ( l + k-^T 1 / 3 |r Jv dm ( 0 

d| 5 |<| a: |/(2TV3) 


h{x,T)= (1 + |£|) 1/2 (1 + \x - £r 1/3 |) n dm{£). 

•7i>m/(2t i /3 ) 

Equation (14.111) now follows from the preceding arguments and the elementary 
estimates 

\h(x,T)\ <c(N) (l + |a;|) -JV 

|/2(*,r)|<c(JV)(l + ^L) ^ . 

This implies (14.81) . 

(4) We claim that the map 

R 2 x R. —> B(L f '( C)) 

(X,T) ^(1-T^r)- 1 

is continuous. If so, it is continuous on the compact set {(x,r) : \x\ < R, |r| < T} 
hence bounded there. By the second resolvent formula and duality, it suffices to 
show that the map 

l 2 xl —> B{L f '(C)) 

( X,T) I y T' x t 

is continuous. But, viewed as multiplication operators on L r (C), e lTS ^'' x ,T —> 
£ itS( ■ ,x,t) ag (x,t) in the strong operator topology. Since 

y' ; ,—T' = (e irS (' ,X ' ,T '^ — e irS (' ’ x,t A W 

for a fixed compact operator W, it follows that 

\\ T x',t' as 0 E '>' r ') -» (x,t). 

This proves the required continuity. □ 

Given the uniform resolvent bounds, we can use the second resolvent formula 
and the continuous dependence of T XjT on its parameters to prove various continuity 
results about the resolvent. 
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Lemma 4.4. Suppose that s £ L 2 (C). For any r > 2, 
(i) The mapping 


I 2 xR —> B{L f ) 

(x,t) h-x (J — T X)T ) _1 


is continuous. 

(ii) The mapping 


is continuous. 


L 2 (C) —5- C°(R 2 x M,£(I7)) 

S l-A ((x, r) l-A (/ - T XiT ) _1 ) 


Proof, (i) This was already proved in the proof of Lemma 14.31 

(ii) Fix ( x,t ) and write T XjT = T XiT (s) to emphasize the dependence on s. From 
the second resolvent formula 

(/-T XjT (s 1 ))- 1 -(/-T ;CiT (s 2 ))- 1 

= (I~ T Xt r( s,))- 1 [T x , r ( Sl ) - T XiT { s 2 )] (/ - T^Si))” 1 

and the fact that 


sup ||T XjT (si) - T Xjr (s 2 )|| B ( L r) < C ? ||si - s 2 || 2 

(i,r)£l 2 xE 


we easily deduce that 


(I~T x ,r( S)Y 


B(L f ) 


is bounded uniformly in (x, t) for s 


in a small metric ball in L 2 whose radius depends on the center but is uniform in 
( x,t ). For Si and s 2 in such a metric ball B we may estimate 


sup 

(x,t)£I 2 x 1 


(J-T x , T (s 1 ))- 1 -(J-T x , r (s 2 ))- 1 


B(Lr) 


<C'(B)||s 1 -s 2 || 2 . 


This gives the claimed continuity. 


□ 


We have already established in Lemma (12.91) that, for s £ Xf r: the derivatives 
d™D x (fj.(x, ■ , t)) exist in 17 (C ), provided that (3m + |a|) < n. We are now ready 
to prove the continuity of several key maps in s. 

Lemma 4.5. Given r £ (1,2) 

(i) For any n > 3m + \a\, the maps 

M7(R 2 xR, 17(C)) 

s ((x,t) d^D"(fJ,{x,;T) -1)) 

are continuous. 

(ii) For any n > 3 m + |a| + 2, the map 

x^ r i-a c^r 2 x R) 

s <9™D“g(x,T) 


is continuous. 
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Proof. For any n > 3 m + |a| we will show the map 

^ r -^C(t 2 xl,L f ( C)) 

ShA ((x,r) S™D“[T x , r (s)(l)]) 

is continuous. The argument is similar to the proof of [19, Lemma 5.4]. From the 
computation 


(4.13) 


d r "rrT x . T i = p k 


e its {-2ik 2 ) ai {2ik 1 ) a2 [i{k 3 + fc 3 )] m s( •) 


and (EH), we see that it suffices to bound || | • •) — s 2 ( • ))|| r by ||si — s 2 || A . e , 

where I = 3 m + |a| < n. This is immediate from the definition. 

Part (i) is the same as m Lemma 5.5] except here we keep track of the continu¬ 
ous dependence on the scattering data. We need to check that the expression (12.91) 
defines a continuous C(M 2 x R; L r (C))-valued function of s. By Lemma [4.414 1). it 
suffices to show that suPi/ (with / given by (12.101) 1 has the same property. A 
typical term of P k f takes the form 

P k [s{-)e its k 0 (i ((-) 3 + (~) 3 )Y {did^i) (x,-,r)] 


where \f3\ + I < m — 1. Assuming inductively that the derivatives d x d^p are 
continuous C(M 2 x R, L r (C))-valued functions of s, we can now use (12.41) with 
q = r, p = r, and s = 2 to obtain the required continuity. Thus it remains to prove 
that p(x, k) — 1 is a continuous C(R 2 x R; L r (C))-valued function of s. This is an 
immediate consequence of Lemma 14. 4f iif and equation (14.131) . 

Part (ii) follows directly from part (i). The proof is the same as in Music m 
Lemma 5.6]. The result follows from equation (11. 91) . part (i), and the fact that 
s G A4 +1>r C L\. □ 


5. Symmetries of the scattering transform 


In this section we prove that q(x) is real if and only if t(fc) = t(— k). This 
symmetry was proved in the negative energy case by Grinevich and Manakov in 
ITO'I where they also sketch an argument for the symmetry in the zero energy case. 

Here we will write this argument directly in the language of the zero-energy 
case. First, we note that Lemma [4711 which allows us to prove relevant formulas for 
smooth potentials and conclude that these formulas continue to hold for q £ L([(R 2 ) 
by continuity. 


Lemma 5.1. If q £ L([(R 2 ) for p £ (1, 2) and p > 2/p' is real and k and —k are 
not exceptional points, then t(fc) = t(— k). 

Proof. First assume that q £ R 2 ) so that the corresponding p(x, k ) is smooth. 

Consider the form 

w = e k d x p(x , —k)dx + e k d x p{x, k)dx 


where 


dx = dx i + idx 2 , dx = dx\ — idx 2 - 


dio = e k ^( d x + ik)d x p{x , k) — ( d x — ik)d x p(x, —k 
= e k (q{x)p(x, k) — q(x)p(x, —k)^j dx A dx. 


dx A dx 


We have 
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Hence, by Stokes’ Theorem, 

t(k) — t(— k) = / ek(x)q(x)/j,(x,k) — ek(x)q(x)n(x,—k) dm(x) = 0. 

Jr 2 

To obtain the result for general q, we approximate q £ L^(R 2 ) by a sequence {q n } 
from C£°(R 2 ) and appeal to Lemma 1-1. 21 □ 

Now, assuming that the scattering data t (k) has this symmetry, we prove that, 
if q(x) is computed using the reconstruction formula (11.51) . then q(x) = q(x). For 
simplicity, we also assume that t(k) =0 in a neighborhood of the origin. We follow 
the proof presented in [9] Theorem 3.4], just rewriting it in terms of the zero energy 
problem. 

Lemma 5.2. If t £ C£°(C) with t(fc) = 0 for |fc| < e for some e > 0 and t(fc) = 
t(— k) then q(x) defined by equation (11.51) is real. Moreover, the solution /i(x, k) of 
the d-problem m satisfies 

d x (d x + ik)ia(x, k) = q(x)p.(x, k ) 

in distribution sense. 


Proof. Consider the real differential form 

n{x,k)n(x,-k) p(x, k)p(x, -k) - 

uj = - ; - dk H-=- dk 

k k 


Using the symmetry t (k) = t (—k) and (11.41) . it is not difficult to see that w is a closed 
form. It now follows by Stokes’ Theorem applied to the region R~ x < |fc| < R , the 
large-fc asymptotic behavior of /.i and / 2 , and the identities 



true for any simple closed contour 7 , that 


(5.1) 


d(x, 0) 2 


p(x, 0 ) 


Thus, fi(x , 0) is either purely real or purely imaginary. Next consider the differential 
operators 


Pi*l> = -d x {d x +ik)ip + qip 
P 2 iP = -d x (d x -ik)ip + qif 

where q is defined by the expansion (12.111) to be q = id x a\. Let xi = Pit 1 an d 
Xi = Pip- We need to show xi = X 2 = 0. From the expansion (12.111) . we have 


lim x'i( a b^) = 0. 

| fc|—>-00 

By (E3D, we have 

(5.2) Xi( x > 0 ) 2 — X2(x, 0 ) 2 = 0, 

and Xi an< i X 2 satisfy 

(dkX 1 ) = e- x (k)s(k)x 2 (x,k) 

{dkX 2 ) {x, k) = e x (k)s(k)xi{x, k). 
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where, for each fixed x, 

lim X 2 (x, k) = lim id x a\ + q + 0 (J/cp 1 ) = id x di(x ) + q{x). 

|fc|—>-oo | k | —yoo 

It is not (yet) clear that id x di{x) + q{x) = 0 but we will prove this using the 
condition (15.21) . 

To this end, consider the one-form 


Xi(k,z)xi(~k,z) dk \ 2 {k. z)\ 2 {-k, z) < - 


V = 


k k 


A computation analogous to the one for u> together with (15.21) shows that rj is a 
closed form and that 

Xi{x,oo ) 2 -X2{x,oof = xi(x,0) 2 - y 2 (a;,0) 2 

where Xi{ x > °°) means linrni^oo Xi( x > k) for i = 1,2. By (15.21) . the right-hand side 
is zero, and also XiC^j 00 ) = 0- We may conclude then that X'Ax- oc) = 0, thus 
q = id x ai =q as desired. □ 

Combining the above result with Lemma 14.51 we have: 

Corollary 5.3. If s £ A’| r (C) forr £ (1, 2) and e > 0 then q(x) defined by equation 
(11.51) is real. Moreover, the solution fj,(x, k) of the d-problem (11.41) satisfies 


d x {d x + ik)n{x, k) = q(x)n(x, k) 


in the sense of distributions. 

It is now simple to prove the interesting fact that all scattering transforms in 
X* r come from critical or subcritical potentials. 

Proposition 5.4. If s £ X.f r then q(x) = [Qt](x) is critical or subcritical. 

Proof. This result follows easily using an approximation argument. Define q f (x) = 
Q(t e )(x) where t e (k) = Xe(k)t(k) and Xe is a smooth radial function which is 1 for 
|fc| > e and 0 for |fc| < e/2. Nachman [2T], Theorem 4.1] proves that the solutions 
H e (x,k) to equation (11.81) satisfy inf^fc \[/, € (x, k)\ > 0 and from (15.11) /z e (x, 0) is 
either real or imaginary. By multiplying this by the appropriate constant and using 
Corollary 15.31 we have cfj, e (k, 0) is a positive solution to the Schrodinger equation 
with potential q e . By Lemma 14.51 q e £ C°(R 2 ), so by [7) Theorem 2.12] q e is 
critical or subcritical. By Lemma 14.51 on any compact subset 11 C l 2 we have 
lmq-jo | q t (x) — g(x)| = 0. By Definition 11.11 a potential q is either critical or 


subcritical if the associated quadratic form (11.101) is positive. That is, for every q e 
and any i/j £ C^°(R 2 ) 


Taking limits preserves the non-negativity and the result follows. 



□ 


The following identities are a consequence of n(x, k, r) solving 

dx(d x + ik)n = qfi 


point-wise as a map k —> L r and the large-fc expansion (12.111) . 
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Corollary 5.5. Assume s £ Xf l+lr an d define q = id x a\. Then, the following 
identities hold: 


(5.3) 


id x u n — O x O x cL n —i ( iO x cLi)ci n —i , 


which for n = 2 simplifies to 

(5.4) id x a 2 = d x 
Additionally we have 

(5.5) id x a 2 = d x 


{^-d x ai + . 

i^-d x ai +iy^ . 


Proof. Identity (15.311 follows by plugging in the large- k expansion (12.111) into (11.81) . 
We get identity (15.51) by applying Liouville’s theorem to the analytic function ia 2 — 

a 2 

d x a\ + i-fi and noting that a\, a 2 , and d x a\ are bounded. □ 


6. Solving the Novikov-Veselov equation 


In this section, we show that q(x, r) defined by (11.91) solves the NV equation 
at zero energy. We will assume s £ X 3r - Note that, for any s of this form, the 
function q(x, t) defined by (11.91) is real by Lemma 15.31 and the fact that the map 
t e lT ( k3 + k )t preserves the symmetry t(fc) = t(—k) and the norms. For such 
transforms, we will prove that q(x, t) does indeed solve the NV equation. It will 
then follow from Lemma 12.61 that the inverse scattering method produces a global 
solution for critical or subcritical initial data q 3 with q 3 £ Wp ,p (K 2 ) for some 
p £ ( 1 , 2 ), and p > 1 . 

This section exploits ideas from Section 3 of the paper of Croke, Mueller, Music, 
Perry, Siltanen, and Stahel 0 , where the Manakov triple representation of the NV 
equation m is used to study the inverse scattering method. In order to show that 
q(x, t) as defined by (11.91) in fact solves the NV equation, we need an equation of 
motion for p{x,k,r). 

We will show directly from the 9fc-equation for p that p{x, k , r) satisfies an 
explicit equation of motion. It then becomes a matter of careful computation to 
show that the time derivative of m satisfies the NV equation. 


Lemma 6.1. If s £ T”! r for r £ (1, 2) and e > 0 then the function p(x, k , r) defined 
by (ED satisfies the equation 




{d T — ikr)p = d x + ( d x + ik ) 3 — 3 u{d x + ik) — 3 ud x p{x , k ) 


where u = id x a\. 


Proof. We have the two equations 

d x (d x + ik)p{x , k) = q(x)p(x, k) 

and 

d k p(x, k, t) = s(fc, t )e- x (k)]J 
Expanding p(x, k, r) up to order fc 3 , we have 

ai(x,r) a 2 (x, r) a 3 (x,r) 


( 6 . 1 ) 


p(x, k, t) = 1 + 


k 2 


k 3 


+ o(k 3 ). 
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By commuting the following specially chosen differential operators through the dk 
equation, we obtain the three identities 

dk(d T - ik 3 )n = s(k, t) e- x (k)(d T - ik 3 )/x, 


TT r-^3 


-pr3 


d k [d x + ( d x + xfc) dM = s(fc,r)e_ x (fc) d x + ( d x + ik) 3 


M) 


and 


dk[— 3ud x — 3u(d x + ik)]/.i = s(fc, r)e_ x (fc) [— 3ud x — 3 u(d x + ik)] /x. 


We can combine these in such a way to get a formula for d T fi. Adding the three 
identities we conclude that 


T = 


/T t3 


d T — ( d x + d x + 3 ikd x — 3 k d x — 3 u(d x + ik) — 3ud x ) 




satisfies 

dk'b = s (k,r)e^ x (k)'If. 

We will prove that dx = 0{k~ x ) so that Theorem 12.31 implies that $ =0. Plugging 
EH into the expression for T, it is easy to see that we get no terms of order k 2 or 
higher. Collecting all terms of order k 1 we find 

—3kd x ai — 3 iku = 0 

which is zero by our choice of u. Collecting terms of order k° gives us 


( 6 . 2 ) 


3id^ai — 3d x a\ — 3iua\ = 3 id x 


d x a i + ia 2 - -a 2 


which is zero by identity (1531) . Thus, the order k° terms in equation EH are zero. 
This finishes the proof that SP = 0. □ 


We are now ready to prove that q(x,r) = id x ai(x,r) solves the Novikov-Veselov 
equation. This follows by looking at the terms of order fc -1 after plugging in the 
asymptotic expansion into the evolution equation for /x and using the identities 
from Corollary 1 5.5 1 

Corollary 6.2. If s £ X§ r then the function q(x,r) satisfies 

d T q = d 3 x q + d x q - 3 d x (uq) - 3 d x (uq). 

Proof. We expand \I> at order /c -1 to get 

d T a\ = d x a\ + d x a\ + 3id^.a2 — 3d x a 3 — 3ud x a\ — 3iua 2 — 3ud x a\. 
Applying the operator id x to both sides we get 

d T q = d x q + d x q - 3d 2 d x a 2 - 3id x d x az - 3d x {u 2 ) + 3 d x (ua 2 ) - 3 d x (uq). 

Now we use equation ESP with n = 2 to get 

d T q = d 3 x q + dlq - 3d^d x a 2 + 3d 2 d x a 2 - 3d x {qa 2 ) - 3d x {u 2 ) + 3d x (ua 2 ) - 3 d x (uq). 

d T q = d 3 x q + d%q - 3 qd x a 2 - 3 d x (u 2 ) + 3 ud x a 2 - 3d x (uq). 

The final step is to show 

—3qd x a 2 + 3ud x a 2 - 3 d x (u 2 ) = -3 d x (uq). 
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We use 


equations (15.411 and (15.511 to get 


3 iqd x (^-d x ai - 3iud x d x a\ + iy- 


-3d x (u 2 ) 

= —3 qd x u + 3 ud x q — 6 ud x u 
= —3qd x u — 3 ud x q 
= -3 d x {uq) 

□ 


The main result of our paper now follows trivially. 

Proof of Theorem \1.2\ By Lemma [2.61 we have s G Xf r when q$ G WA p (R 2 ) for 
p G (1,2), p > 1, e > 0, and r G (p',oo). By Corollary 16.21 q(x,r) solves the 
Novikov-Veselov equation. By Lemma T4.5I part (iii) q(-,r) G (7(R 2 x R) and hence 
q{x,r) q(x,0), where g(ai,0), the inverse scattering transform of s(fc, 0), is the 
initial datum. □ 


7. Conclusion 

We have shown that we can solve the Novikov-Veselov equation via inverse scat¬ 
tering for a wide range of potentials. The most significant problem left in zero- 
energy inverse scattering is the handling of supercritical initial data, for which the 
exceptional set need not be empty and the scattering transform may have non- 
integrable singularities of codimension one (see, for example [201 ). In addition, 
Taimanov and Tsarev [26, §4] have constructed supercritical initial data for the 
NV equation corresponding to a solution that blows up in finite time. This solu¬ 
tion currently lacks in interpretation in terms of inverse scattering, and such an 
interpretation would be of considerable interest. 

Another open problem concerns the decay of solutions with subcritical initial 
data. Lassas, Mueller, and Siltanen m prove that certain reconstructed critical 
potentials have the decay |(Qt)(x)| < C(x)~ 2 . We would like to know that this 
result, or a similar one, holds for subcritical potentials as well, but the singularity 
at k = 0 in the scattering transform poses additional difficulties. 
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